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I.  Introduction  and  Review 

A.  Relation  to  other  computer  simulation  projects 

There  have  been  several  past  and  current  computer  simulations  of  stripe 
damage  caused  by  electromigration,  some  of  which  were  listed  in  the  references 
accompanying  my  original  proposal.  Most  of  these  programs  have  been  run  many 
times  to  afford  statistical  information  on  the  half-life  to  failure  and  the 
corresponding  random  spread  parameter.  Thereby  they  have  been  able  to  afford 
considerable  valuable  information  on  the  influence  of  grain  size  and  of  stripe 
dimension  and,  in  an  approximate  way,  to  take  into  account  the  effects  caused 
by  current  and  temperature  inhomogeneities . 

Our  project  differs  from  these  other  simulations  in  two  respects:  (1)  we 
have  tried  to  make  our  model  somewhat  more  realistic  and  detailed  and  (2),  as 
a  consequence  of  the  increased  complexity,  we  shall  probably  be  limited  to 
only  those  initial  stages  of  damage  before  the  effects  of  current  and 
temperature  inhomogeneities  become  important.  Although  well  developed 
statistical  results  will  be  difficult  to  come  by,  we  can  explore  the  influence 
of  texture,  annealing  and  surface  condition  on  damage  by  electromigration. 

The  approach  should  provide  a  useful  alternative  to  other  computer 
simulations . 

B.  Review  of  prior  results  in  this  area: 

The  initial  work  in  setting  up  this  project  was  concerned  with  developing 
a  suitable  two-dimensional  network  to  simulate  the  {111}  texture  of  the 
metallized  stripe.  For  this  a  Voronot  procedure  was  applied  and  the  grains  of 
the  structure  were  assigned  random  oi ientations .  The  grain  boundary  behavior 
followed  the  dislocation  model  for  misorientation  angles  less  than  32°  and  was 
assumed  to  be  isotropic  for  greater  misorientations .  On  the  basis  of  this 
model  expressions  were  developed  for  both  the  grain  boundary  dif fusivities  and 
the  grain  boundary  tensions. 
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The  next  step  was  to  relax  the  grain  boundary  tensions  by  displacing  the 
network  vertices  appropriately.  In  Figs.  1  and  3  we  show  Voronoi  pattern  for  a 
wide  stripe  A  (15  pm)  and  a  medium  stripe  B  (10  pm).  In  adjoining  Figs.  2  and 
4  are  shown  the  same  stripes  after  equilibration  of  g.b.  forces.  This  process 
could  be  thought  of  as  a  form  of  annealing.  It  has  bearing  on  the 
electromigration  damage  process  in  that  the  condition  for  stress  relaxation  by 
vertex  motion  along  the  axis  parallel  to  the  current  is  closely  related  to  the 
equation  for  mass  accumulation  (or  deficit)  at  the  vertex  in  question. 

In  proceeding  now  to  the  actual  electromigration  one  writes  out  the 
diffusion  equation  appropriate  to  conditions  of  forced  motion  from 
electromigration.  Working  toward  a  solution  of  this  equation  leads  to  two 
approximations : 

1.  Because  the  time  constant  of  the  network  segments  (L^/D)  were  so 
short,  it  was  reasonable  to  bypass  the  transient  solutions  and  move 
directly  to  steady-state,  or  quasi  steady-state,  solutions. 

2.  It  became  apparent  that  the  principal  effect  of  the  electromigration 
on  the  mass  distribution  in  the  grain  boundaries  came  about  directly 
at  the  vertices  and  its  effect  secondarily  on  the  concentration 
distributions  in  the  segments  was  small.  Accordingly  the  electro¬ 
migration  term  in  the  diffusion  equation  was  dropped  and  only 
considered  right  at  the  vertices. 

At  each  vertex  one  evaluates  Fv  =  EgDgAo  cos<f>0  Mq  where  the  summation  is 
over  the  various  segments  that  meet  at  the  vertex,  D0  is  the  diffusion 
coefficient  for  each  segment  as  determined  by  the  mlsor ientat ion  between  the 
grains  on  either  side  of  the  grain  boundary,  and  the  angles  <f>0  are  between  the 
direction  of  the  grain  boundary  and  the  electric  current.  The  number  of  atoms 
per  unit  length  of  grain  boundary  is  MQ .  The  constant  A0  is  |e|Z*E/kT  where 
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e  is  the  electron  charge,  E  the  electric  field  and  T  the  temperature.  Z*  is 
the  effective  charge  number  for  electromigration  which  is  around  -20  for 
aluminum.  If  Fv<0,  the  vertex  is  a  deficit  vertex  (-)  and  conversely  if  Fv>0, 
the  vertex  is  an  excess  vertex  (+). 

The  (+)  or  (-)  vertices  are  treated  quite  differently.  The  net  deficit 
at  the  (-)  vertices  is  considered  to  contribute  directly  to  a  growing, 
circular  cylindrical  void.  At  the  (+)  vertices  the  matter  can  not  pile  up  at 
a  point  but  must  diffuse  back  along  the  grain  boundaries.  As  a  result  the 
( — )  segments  maintain  the  uniform  grain  boundary  concentration  according  to 
our  model.  For  the  (+-)  segments  the  concentration  goes  to  zero  at  the  (-) 
end  but  may  build  up  at  the  (+)  end.  There  is  build-up  at  both  ends  of  the 
(++)  segments. 

We  have  been  concerned  with  electromigration  under  two  extreme  surface 
conditions:  free  surface  and  the  surface  covered  with  a  passivating  layer 
which  prevents  flow  up  out  the  grain  boundary.  These  two  situations  will  be 
discussed  respectively  in  the  next  two  sections. 

II.  Electromigration  under  the  Free  Surface 

According  to  the  model  the  excess  material  from  the  (+)  vertices  diffuses 
down  the  grain  boundaries  'and  up  out  of  the  grain  boundaries  to  cover  the 
surface  quite  smoothly.  This  situation  can  be  expressed  mathematically  simply 
by  solutions  which  are  steady  state  within  the  grain  boundaries.  The 
solutions  must  satisfy  the  Laplacian  and  will  be  of  the  form  cosh  (or  sinh) 

1TZ 

yx  cos  where  x  is  measured  along  the  g.b.  and  z  is  measured  from  the 

substrate  to  the  free  surface.  Accordingly  the  surface  deposits  centered  at 
the  (+)  vertices  grow  with  time  as  do  the  voids  centered  at  the  (-)  vertices. 

We  summarize  the  mathematical  treatment  next,  starting  with  the 
differential  equation  for  mass  movement, 
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where  Ma(x,z,t)  is  the  amount  of  material  in  the  o-graln  boundary  per  unit 
length  (x)  of  boundary.  The  z  is  measured  perpendicular  to  the  surface  from 
the  substrate  to  che  free  surface.  The  subscript  o  denotes  the  segment  under 
consideration.  D0j_  means  the  diffusion  coefficient  perpendicular  to  the 
dislocations  in  the  grain  boundary,  i.e.  along  the  surface,  and  D a j|  is  for 
motion  parallel  to  the  dislocations,  i.e.  at  right  angles  to  the  substrate. 

Let  us  write 

Mo(x,z,t)  =  Mq  +  m(x,z,t)  (-2) 

where  Mq  is  the  equilibrium  number  of  atoms  per  unit  length  of  grain  boundary, 
or  N6d  where  N  is  the  number  of  atoms  per  unit  volume,  d  is  the  thickness  of 
the  stripe  and  6  is  the  width  of  the  grain  boundary  or  roughly  one  atom  dimen¬ 
sion.  The  grain  boundary  segments  can  be  categorized  by  the  nature  of  their 
end  points,  ( — ),  (+-)  and  (++) •  Because  of  void  formation  there  is  no  back 
diffusion  long  the  ( — )  segments  and  m0  =  0.  For  the  other  segments  we  are 
interested  in  the  time  independent  situation  after  the  transient, 


m(x,z,“>)  =  aQ(x)  ~  cos(  irz/2d)  (3) 

The  form  of  the  z-dependence  is  chosen  to  give  all  matter  flow  at  the  surface 

and  none  at  the  substrate.  The  equation  for  m(x)  now  becomes 

ra"(x)  -  A  m'(x)-Cm(x)  *  0  (4) 

a 

YjX  Y2(x-L) 

with  solutions  m(x)  ■  a  e  +  a  e  where  (5) 
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A  A  2  ~  0  „  2 

Y  1,2  -r  1  «2~>  +  c)  and  c  =  d^-<2T-) 

*  ol 


Since  D|j>D^  in  every  case  and  d~ly  while  A~  O.iy,  A/2  is  usually  much 
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2 

smaller  chan  C  so  that  y^  «  -y^  =  ~Y*  This  is  consistent  with  our  second 
approximation  as  stated  earlier. 

We  concentrate  now  at  one  vertex  mv  and  consider  in  turn  each  connecting 
segment  as  indexed  by  a.  For  each  segment  there  are  two  equations,  taking  x  = 
0  at  the  vertex  in  question. 
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where  mva  is  the  value  of  m  at  the  other  end  of  the  o-segment.  For  the  (+-) 

— yl  — 2 vL 

segments  m^  a  0,  e  and  ra^  *  a^Cl-e  ). 

The  magnitudes  of  the  mv,  or  a's  are  determined  by  the  condition  at  the 
(+)  vertices  that  the  diffusion  flux  carries  away  the  excess  material 
deposited  by  electromigration 
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Based  on  eq.  (5). 
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Hence 
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We  now  solve  for  mv 

m  =  [F  (1  -  e“2yL)  +  Z  2e  0  0  m  D  ,y 
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To  solve  equations  12  there  is  the  initial  difficulty  that  the  mv a  are 
unknown.  This  can  be  overcome  by  iteration  on  the  computer.  Let  the  index  n 
designate  the  ra  value  at  the  end  of  the  nth  stage  of  iteration. 
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We  begin  the  first  step  with  m  =0. 

voo 

In  the  course  of  each  step  in  the  iteration  process  there  will  be 

occasion  to  note  for  the  (+-)  segments  a  rate  of  deposit  at  the  (-)  ends  of 

_2 

2DomoYa  exP(“YaL,a)/(l-e  Lyo)  whlch  tends  to  compensate  the  deficit  caused  by 
the  original  Fv  <  0.  We  represent  the  missing  material  by  circular  cylinder. 
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the  quantity  inside  the  {  }  turns  positive,  the  vertex  changes  from  (-) 

(+)  and  must  be  so  considered  +  in  subsequent  iterations  with  Fv’  =  {  }. 

The  result  of  this  electromigration  under  a  free  surface  is  illustrated 
Figs.  5-8. 

The  technique  we  have  employed  is  to  represent  the  deficit  cylinders  by 
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open  ci.rcl.es  scaled  Co  size  and  the  build-up  regions  by  cross  hatchings  whose 
amplitudes  are  in  proportion  to  the  excess  material  in  the  grain  boundaries. 
There  was  a  need  to  choose  an  index  to  measure  the  extent  of  damage  for  each 
stripe.  For  this  we  took  fc,  the  fraction  of  stripe  volume  converted  to  void 
by  the  electromigration  action.  The  fc  -  values  for  Stripes  A.  &  B  before  and 
after  equilibration  are  given  in  Table  1.  Note  the  extent  to  which 
equilibration  has  improved  stripe  resistance  to  electromigration  damage. 

III.  Electromigration  Under  a  Passivating  Layer 

This  situation  is  just  the  opposite  extreme  from  that  of  the  free 
surface  treated  in  the  preceding  section.  Any  escape  of  the  material  in  the 
grain  boundary  to  the  surface  is  prevented  by  the  passivation.  The  pressures 
in  the  grain  boundaries  connected  to  excess  vertices  thereby  increase  to 
limiting  values,  assuming  no  rupture  of  the  passivation.  This  eventually 
attained  steady  state  occurs  when  the  electromigration  influx  at  any  excess 
vertex  is  balanced  by  the  outward  flow  to  nearby  (weaker)  excess  vertices  or 
to  deficit  vertices. 

The  build-up  period  described  above  has  proved  very  difficult  to  model 
satisfactorily.  We  have  made  several  attempts  which  have  failed.  Currently 
we  have  developed  an  approach  which  looks  promising,  but  testing  with  a 
suitable  program  has  been  delayed  as  a  new  programmer  is  learning  the  ropes. 

This  approach  is  built  around  the  use  of  a  particular  function  which  is  a 
solution  of  the  diffusion  equation  in  the  segments. 

C(x, t)  -  (x/L)  +  6D  t  =  (x/L)  +  6  —  (t/xr)  (11,1) 

L 

where  L,  D  and  t  are  respectively  the  segment  length,  the  diffusity  in  the 
segment  and  the  segment  time  constant  t0  =  Do/L0^.  Accordingly  we  start  the 
build-up  timing  after  a  short  transient  period  and  assume  that  the 
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m-distr ibuticm  in  the  (+  -)  segments  is  given  by 

ma(x* t)  =  Wo^c)  co(x»t)  (II, 2a) 

and  in  the  (++)  segments  by 

m0  (x,t)  =  pL(t)  Ca{(L-x),  t}  +  pR  Ca  (x,t)  (II, 2b) 

where  subscripts  L  and  R  refer  to  the  right  and  left  ends  of  the  segment 
and  the  p  are  presumably  slowly  varying  functions  of  time.  The  presence  of 
the  p's  does  of  course  prevent  the  m  functions  from  being  solutions  of  the 
diffusion  equation^  particularly  as  the  function  of  the  p  is  to  bring  the  id's 
to  steady  state  values  and  to  stop  the  linear  increase  with  time.  We  believe, 
however,  that  this  is  not  a  serious  weakness,  since  one  can  devise  a  solution 
which  is  a  series  of  functions,  dominated  by  a  first  term  which  saturates.  It 
remains  now  to  invoke  the  condition  established  for  the  Fv. 

+  -  D  ++  D 

F  =  l  (3  +  6  — )  p  (t)  +  )  [t^-  [3  +  —  (p  -  p  J]  (11,3) 

v  L  L  t  '  v  L  1 L  1  pv  T  V  vcrJ  v  ’  ' 

a  a  a  a  a  a 


At  the  end  of  the  transient  (t  =  o) 

Uv(0)  =  mv  (O)  =  F v/l  3  (D/Lo  (11,4) 

a 

After  a  short  interval,  At, 


p  (At)  =  {F  +  16  {D/L  (At/t  )  p  }/l  (D/L)  (3  +  IL^I) 

v  v  u  o  ovoJ  a  x 

o  o 

Here  the  pv0  are  the  p-values  on  the  (+)  vertices  which  are  linked 

For  them  we  approximate  with  pv0  (o)  =  my(o)  initially.  It  may  be 

foLLow  later  with  one  stage  of  iteration,  using  Eq .  II,  5  a  second 

the  pvg(^(At)  on  the  right  hand  side  to  get  improved  pv(At).  For 

time  steps  it  shouLd  be  sufficient  to  extrapolate 


(11,5) 

to  vertex  v. 
advisable  to 
time  with 
succeeding 


[n+-H  ,/Tr  n  n 

p'  2J  =  1/2  [  p  -p  +  p 

VO  1  VO  VO  VO 


to  be  used  in  the  extension  of  (II,  5)  to  long  times. 


(11,6) 
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Obtaining  my,  now  that  the  Uy  are  determined,  runs  into  a  difficulty  with 
the  C-functlons  as  introduced  in  (It,  2).  Our  method  insures  that  the  ^  is  a 
function  of  v  only  but  this  is  not  true  of  C0  where  one  sees  a  a  dependence  in 

3 ' 

1  11'  1 

the  (t/xa).  To  bypass  this  difficulty  we  replace  —  by  —  =»  —  jf  — . 

a  v  a  a 

With  the  use  of  this  crude  approximation  we  can  now  write 

n>v  (t)  -  uv(t)  U  +  6  t/rv}  (11,7) 

Ac  each  step  in  the  time  progression  one  must  also  calculate  the  niy  from 
II,  2  and  also  the  spill-over  at  the  (-)  vertices.  For  these  vertices  the 
void  volume,  fty,  is  given  by 
n  +  — 

W  -  l  I-  i  6<V/V>  -  Fv)  <4t>  -  (11,8) 

n  =  l  a 


n 

where  t^  *  \  (At)^-  and  F  <  0  for  (-)  vertices 

The  Ciisc  term  in  the  curly  brackets  comes  from  the  spill-over  effect 
(S-0).  If  the  S-0  term  becomes  bigger  than  -  Fv  the  void  may  fill  up.  If 
this  happens,  the  program  should  recognize  that  the  vertex  has  changed  from 
(-)  to  (+)  with  an  effective  Fv "  =  -  {  }  of  equation  (II  7).  It  will  be 

important  to  get  some  measure  of  the  extent  to  which  the  original  number  of 
voids  is  reduced  by  this  filling  process. 

What  form  will  the  riiy  take  at  long  time? 

One  can  obtain  from  Eq .  II  5  a  set  of  linear  relations  for 
mv  =  c/ T  as  defined  by  Eq .  II  7 

T  +  ao 

These  relations  are 


++ 

=ny  (-)  l  (D/L)o  -  Fv  +  l  (D/L)  m  ( -) 
a  a 


(11,9) 
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by  transforming  Eq.  (II  5).  This  shows  that  it  is  possible  to  obtain  the 
final  values  for  mv  without  following  through  the  detailed  development  in 
time.  A  method  of  checking  the  detailed  treatment  Is  thereby  provided. 

This  procedure  appears  promising  in  that  it  seems  to  have  avoided  the 
singularities  and  divergencies  which  wrecked  earlier  attempts  at  developing  a 
treatment  of  the  passivated  layer  case.  Actual  implementation  is,  however, 
held  up  temporarily  while  a  new  programmer  is  learning  the  operation  of  the 
project . 

IV.  Future  Plans  (from  Jan.  1988  report) 

1.  If  the  procedure  outlined  in  the  previous  section  works  well  for  the 
strongly  passivated  stripe,  we  should  learn  how  long  it  takes  to  reach  the 
steady  state  and  at  what  build-up  levels,  what  fraction  of  the  voids  become 
filled  and  what  the  influence  of  stripe  dimensions  and  stripe  annealing  might 
be.  For  a  truly  realistic  picture  one  needs  to  make  provision  for  the  rupture 
of  the  passivity  layer  and  the  formation  of  a  hillock.  Presumably  this  will 
occur  at  a  predetermined  stress  level  or  perhaps  value  of  mv.  From  there  on 
out  rav  at  this  vertex  Is  zero.  The  Fv  from  electromigration  goes  directly  to 
increasing  the  "whisker"  and  has  no  effect  on,  m(x).  Likewise  now  the 
diffusional  flux  is  inward  at  the  vertex  in  question  and  the  general  tendency 
for  void  formation  in  the  neighborhood  is  increased,  as  has  been  observed 
experimentally.  This  occurs  mainly  because  of  the  reduction  in  "spill-over". 

2.  We  also  want  to  explore  the  effect  of  surface  electromigration  which 
will  cause  void  motion  in  the  opposite  direction  from  the  mass  motion.  For 
simplicity  we  take  the  void  cylinders  to  be  initially  circular  and  assume 
that  surface  diffusion  maintains  this  shape  as  that  of  minimum  surface  energy. 
The  mass  transport  from  one  half  of  the  void  surface  to  the  other  is  given  by 
ZDgAjjd  as  calculated  at  those  surface  elements  of  the  void  parallel  to  the 
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current.  Here  Dg  is  the  surface  diffusion  coefficient.  Dividing  the  transport 
by  the  number  of  atoms  missing  in  the  void  gives  the  void  velocity,  vc . 

2 

V  =»  2D  A  /NirR 
c  so 

Apparently  the  smaller  voids  move  more  rapidly  and  so  tend  to  catch  up  and 
coalesce  with  the  larger  voids.  The  process  of  surface  electromigration 
also  tends  to  distort  the  cylinders  into  cracks  along  the  grain  boundaries,  a 
feature  that  may  be  difficult  to  simulate.  Nevertheless  it  may  be  quite 
important  in  the  deterioration  process. 

3.  Another  feature  which  we  plan  to  consider  is  the  possibility  of  invoking 
grain  growth  more  directly  in  the  annealing  process  than  results  from  the  • 
stress  equilibration  process.  This  would  be  done  by  looking  for  those  grains 
which  may  disappear  altogether  as  a  result  of  surface  tension. 

4.  One  interesting  ramification  of  the  above  might  arise  if  one  postulates  as 
a  result  more  rapid  grain  growth  on  the  surface  of  the  stripe  than  at  the 
stripe-substrate  interface.  As  a  consequence  not  all  grain  boundaries  will  be 
perpendicular  to  the  stripe  surfaces.  Then  some  of  the  voids  may  not  be  right 
circular  cylinders  but  develop  "end"  sections  across  which  rapid  surface 
electromigration  can  occur.  This  effect  could  contribute  very  substantially  to 
void  motion,  but  would  involve  3-dimensional  consideration  which  would 
introduce  formidable  complications. 

5.  On  the  other  hand  the  treatment  of  peaked  electromigration  (or  even  ac 
operation)  appears  straightforward  and  perhaps  well  worth  pursuing. 

The  numerous  effects  and  responses  of  a  stripe  undergoing  electromigration 
are  so  dramatic  and  varied  than  it  is  unrealistic  to  hope  to  incorporate  them 
all  in  a  single  simulation.  We  have  tried  in  this  section  to  propose  only  a  few 
aspects  of  the  general  problem  which  evolve  from  the  particular  model  we  have 
devised  and  are  therefore  susceptible  to  investigation  in  future  studies. 


Table  1 

Values  of  f,. 

Original 

After  Equilibration 

Ratio 

Stripe  A 

3.70% 

2.07% 

.56 

Stripe  B 

3.27% 

2.16% 

.66 
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List  of  Figure  Captions 


1.  Wide  stripe  A  -  Voronoi  construction 

2.  Wide  stripe  A  after  equilibration 

3.  Stripe  B  -  Voronoi  construction 

4.  Stripe  B  after  equilibration 

5.  Effect  of  electromigration  on  Stripe  A  -  before  equilibration 

6.  Effect  of  electroraigration  pn  Stripe  A  -  after  equilibration 

7.  Electromigration  on  Stripe  B  -  before  equilibration 

8.  Electromigration  on  Stripe  B  -  after  equilibration 
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Figure  2  -  Wide  stripe  A  after  equilibration 


Figure  4  -  Stripe  B  after  equilibration 
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